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The topic addressed in this paper was initially motivated by questions resulting from recent lab- 
oratory experiments designed to empirically test and understand advection-dispersion in the presence 
of sharp interfaces; e.g., experiments by [6], [5] [3]. Such laboratory experiments have been rather 
sophisticated in the use of layers of sands and/or glass beads of different granularities and modern 
measurement technology. As a result they have uncovered a convincing empirical foundation for some 
interesting and unexpected phenomena that had escaped prior theoretical notice and explanation. To 
this end it is natural to inquire about the effect of an interface on the stochastic particle motion of 
immersed solutes. From a general mathematical point of view an interface is defined by a hypersur- 
face across which the dispersion coefficient is discontinuous. As is well-known for the case of dilute 
suspensions in a homogeneous medium (e.g., water), perhaps flowing at a rate v, the particle motion 
is that of a Brownian motion with a constant diffusion coefficient D > and drift v. For simple 
one-dimensional flow across an interface, a localized point interface results in a skewness effect that 
explains much of the empirically observed results noted above; see [9], [TO], PQ, [2], [TT] . 



Problem As an illustration of empirical findings, suppose that a dilute solute is injected at a point 
L units to the left of an interface at the origin and retrieved at a point L units to the right of the 
interface. Let D~ denote the (constant) dispersion coefficient to the left of the origin and D + that 
to the right, with say D~ < D + . Conversely, suppose the solute is injected at a point L units to the 
right of the interface and retrieved at a point L units to the left. Which will be retrieved first ? 

For positive parameters D + , D~ , consider a piecewise constant dispersion coefficient with inter- 
face at x = given by 

D(x) = D~T (-00,0) 0*0 + D + l[ 0iOO )(x), x G R. 

Theorem 1 Let D + ,D~ be arbitrary positive numbers, with say D~ < D + . Define = 
s{B[ a '),t > 0, where £?( a *) is skew Brownian motion with transmission parameter a* = D +^ D - , 
and s(x) = V D + xl[ 0tOO )(x) + V D~ xlr^™ (x), x G R. Let T y = inf{t > : Y"/ Q ^ = y}. Then, 

(a) For smooth initial data cq, c(t,y) = E y co(Y t ), t > 0, solves 



dt 2 dy dy ' dy dy 

(b) For y > 0, > t) < ^P^^ > t) < P y (T^ y > t), t> 0. 

Remark This basic result was obtained JE/ in terms of first passage times, however the factor V D~ / \f~D^ 
was not included in the statement of the result there. Related phenomena and results on dispersion 
in this context are also given in flUtf . Jiy, \1 If . In addition, a formula for the first passage time 
distribution for skew Brownian motion was recently obtained in J3]/. In principle, the identification 
of stochastic particle motions can have computational advantages. Results pertaining to Monte- Carlo 
simulations of skew diffusions are described in |W and references therein. 

As illustrated by the examples below, the role of interfacial phenomena is of much broader inter- 
est than suggested by advection-dispersion experiments. However the specific nature of the interface 
can vary, depending on the specific phenomena. We briefly describe three distinct classes of examples 
of phenomena from the biological/ecological sciences in which interfaces naturally occur. 

Example 1 (Coastal Upwelling and Fisheries) Up-wellings, the movement of deep nutrient rich 
waters to the sun-lit ocean surface, occur in roughly one percent of the ocean but are responsible for 
nearly fifty-percent of the worlds fishing industry. The up-welling along the Malvinas current that 
occurs off of the coast of Argentina is unusual in that it is the result of a very sharp break in the 
shelf, rather than being driven by winds. The equation for the free surface r\ as a function of spatial 
variables (x, y) is of the form 

drj r f dh\~ x d 2 r\ 



dy f \dx J dx 2 ' 

where r > and / < in the southern hemisphere, and h{x) is the depth of the ocean at a distance x 
from the shore. In particular, the sharp break in the shelf makes h'(x) a piecewise constant function 
with positive values H + ,H~. The location of the interface coincides with the distance to the shelf- 
break. If the spatial variable y > is viewed as a "time" parameter, then this is a skew-diffusion 
equation, however the physics imply continuity of the derivatives drj/dy at the interface; see [8] and 
references therein. 



Example 2 (Fender's Blue Butterfly) The Fenders Blue is an endangered species of butterfly 
found in the pacific northwestern United States. The primary habitat patch is Kinkaid's Lupin flower. 
Quoting j!4j . "Given past research on the Fender's blue, and the potential to investigate response to 
patch boundaries, we ask two central questions. First, how do organisms respond to habitat edges? 
Second, what are the implications of this behavior for residence times?" Sufficiently long residence 
(occupation) times in Lupin patches are required for pollination, eggs, larvae and ultimate sustainabil- 
ity of the population. Empirical evidence points to a skewness in random walk models for butterfly 
movement at the path boundaries. 

Example 3 (Sustainability on a River Network) The movement of larvae in a river system 
is often modeled by advective-dispersion equations in which the rates are determined by hydro- 
logic/geomorphologic relationships in the form of the so-called Horton laws. In general river networks 
are modeled as directed binary tree graphs and each junction may be viewed as an interface. Con- 
servation of mass leads to continuity of flux of larvae across each stream junction as the appropriate 
interface condition. Problems on sustainability in this context are generally formulated in terms of 
network size and characteristics relative to the production of larvae sufficient to prevent permanent 
downstream removal at low population sizes; see [12] for recent results in the case of a river network. 

The following theorem provides a useful summary of the interplay between diffusion coefficients 
and broader classes of possible interfacial conditions illustrated by these examples. The proof follows 
by a straightforward application of the Ito-Tanaka formula. 

Theorem 2 Let D + ,D~ be arbitrary positive numbers and let < a, A < 1. Define y[ = 
s(B^ a )),t > 0, where is skew Brownian motion with transmission parameter a and s(x) = 

\/S+2 ; l[o,oo) (x) + \/ZF a;l(_ OO)0 ](a;), x G R. Then 

M t = f{Y t {a) ) -\f D{Y^)f"{Y u )du, t > 0, 
^ Jo 

is a martingale for all / G V x = {/ G C 2 (R\{0}) n C(R) : A/'(0+) = (1 - A)/'(0~)} if and only if 

a = a A) = — ; ; . 

A^/lF + (1 - X)VJj+ 



Remark This theorem is a generalization of the results obtained by and JI]/ for the case of advection- 



dispersion problems across an interface described at the outset, where the parameter A = -qt+jp 



and 



a 



n5+- 



'D- 



Definition With the choice of a* = a* (A) given by Theorem 2, we refer to the process Y a * as the 
physical diffusion corresponding to the dispersion coefficients D + ,D~ and interface parameter A. 

Observe that in the application to the coastal up-welling problem one obtains 

a ■ 



f D+ + VD- 



The physical diffusion for this example may be checked to coincide with the Stoock-Varadahn mar- 
tingale in this case; see [15J for the definition of the corresponding martingale problem. Note that 
the answer to the first passage time problem will be exactly opposite to that obtained for advection- 
dispersion experiments under this model. 



We conclude with a result to show that the issue raised in Example 2 relating interfacial condi- 
tions to residence times is indeed a sensitive problem. 

Theorem 3 Let Y a * denote the physical diffusion for the dispersion coefficients D + , D~ and interface 
parameter A. Define modified occupation time processes by 

f + (t) = f l[Y^ > 0}ds, t > 0. 
Jo 

Similarly let f~(t) = t - f+(t), t > 0. Then, 



r + (t) > r - (t) vt > o ^a> 



with equality when A 



fD+ 



Proof Let A(x) = 2Al[ 0iOO )(x) + 2(1 — A)l(_ oo )(x), and define = D(x) / \{x). Consider the time 
change T p (t,oj) defined by 

r-r P (t) i 

/ -7TTT = t, t>0. 
Jo P(B S ) 

Define S p by B(t,S p (u)) = B(r p (t,uj),oj) = Z(t,u). Then the process Z is a diffusion with zero drift 
and diffusion coefficient with interface parameter A = 1/2. Now observe that 



rf(a.)(*) == / o *lK ( ^>0]d<F^)> s 

= f f l[Z s >0}4X 2 d< Z > s 
Jo 

= j\[B(r p (s)>0]^d S 



= 4A 2 r+(r p (t)). 

Similarly Ty {a9) (t) = 4(1 - \) 2 T^(T p (t)). Thus (1 - A) 2 r+ (a ,)(t) = A 2 r- (a .)(t). Now observe that 

/•* 4A 2 

r+ „.,(*) = y o i[z, > o]d s = _r+(r p (t)), 

and similarly for r+^.^t), to arrive at 



(1 - A) 2 w A 2 

The assertion now follows. 



It is interesting to note that under the mass conservation interface parameter A = D + /(D + + 
D~), the particle will reside longer in the region with the faster dispersion rate. While this is to be 
expected for physical experiments of dispersion in porous media of the type described above, it shows 
that the conservative interface condition (defined by this choice of A) is likely not appropriate for 
models of animal movement ! 



Remark A related phenomena in terms of a "modified local time" is described in f^j. The modification, 
denoted with the ~ } refers to an integration with respect to Lebesgue measure in place of quadratic 
variation in the usual mathematical definition of local time and quadratic variation; e.g., see JJ3]- In- 
general, the treatment of dispersion in the presence of interfaces suggests that the physical/biological 
theories are, to the extent possible, naturally based on a modification of local times and occupation 
times in which integration with respect to quadratic variation is replaced by integration with respect to 
Lebesgue measure. In fact it is shown in [2] that this naturally leads to a stochastic determination of 
the physical transmission parameter a* in terms of a continuity condition on the modified local time of 
the stochastic particle motion. This is a probabilistic condition at the particle scale that may be viewed 
as an alternative to the usual macro-scale pde condition of continuity of flux in particle concentrations. 
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ABSTRACT 

Advection and dispersion in highly heterogeneous environments involving interfacial discontinu- 
ities in the corresponding drift and dispersion rates are described through disparate examples from the 
physical and biological sciences. A mathematical framework is formulated to address specific empirical 
phenomena involving first passage time and occupation time functionals observed in relation to the 
interfacial parameters. 



